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ABSTRACT Horizontal gene transfer (HGT) is a major factor in the evolution of prokaryotes. An intriguing KEYWORDS 
question is whether HGT is maintained during evolution of prokaryotes owing to its adaptive value or is environmental 
a byproduct of selection driven by other factors such as consumption of extracellular DNA (eDNA) as DNA 
a nutrient. One hypothesis posits that HGT can restore genes inactivated by mutations and thereby prevent evolution of 
stochastic, irreversible deterioration of genomes in finite populations known as Muller's ratchet. To examine transformation 
this hypothesis, we developed a population genetic model of prokaryotes undergoing HGT via homolo- competence 
gous recombination. Analysis of this model indicates that HGT can prevent the operation of Muller's ratchet structured 
even when the source of transferred genes is eDNA that comes from dead cells and on average carries more population 
deleterious mutations than the DNA of recipient live cells. Moreover, if HGT is sufficiently frequent and soil bacteria 
eDNA diffusion sufficiently rapid, a subdivided population is shown to be more resistant to Muller's ratchet 
than an undivided population of an equal overall size. Thus, to maintain genomic information in the face of 
Muller's ratchet, it is more advantageous to partition individuals into multiple subpopulations and let them 
"cross-reference" each other's genetic information through HGT than to collect all individuals in one 
population and thereby maximize the efficacy of natural selection. Taken together, the results suggest that 
HGT could be an important condition for the long-term maintenance of genomic information in prokaryotes 
through the prevention of Muller's ratchet. 



Evolutionary advantages of recombination have been extensively con- 
sidered, primarily in relation to the evolution of sex in eukaryotes 
(Kondrashov 1993; Maynard Smith 1978, 1998). However, advances 
in genome sequence analyses have revealed that many prokaryotes, 
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whose mode of reproduction is asexual (i.e., no meiosis or syngamy), 
frequently undergo genetic recombination, the process known as hor- 
izontal gene transfer (HGT; Feil et al 2000, 2001; Ochman et al 2000; 
Koonin et al 2001; Gogarten et al 2002; Vos and Didelot 2009; 
Yahara et al 2012). Numerous studies indicate that HGT has been 
instrumental in producing functionally consequential genomic 
changes throughout the evolution of prokaryotes (Lerat et al 2005; 
Dagan et al 2008; Treangen and Rocha 2011). For example, HGT is 
implicated in the evolution of various functional systems and specific 
adaptations such as various metabolic pathways (Pal et al 2005; 
Maslov et al 2009), oxygenic photosynthesis (Mulkidjanian et al 
2006), thermal resistance (Aravind et al 1998; Brochier-Armanet 
and Forterre 2007), antibiotic resistance (Blahna et al 2006), patho- 
genicity (Kelly et al 2009), and others. Moreover, it has been shown 
that HGT can substantially impact microbial community develop- 
ment (Burke et al 2011) and ecological diversification (Shapiro 
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et al. 2012). The prevalence and significance of recombination in 
prokaryotes as well as in eukaryotes point to the crucial importance 
of recombination for the evolution of life in general (Komatsu 1980). 

HGT in prokaryotes occurs though three basic mechanisms: 
conjugation, transduction, and transformation (Levin 1988; Thomas 
and Nielsen 2005; Johnsborg et al. 2007). Conjugation and transduc- 
tion involve transfer of prokaryotic DNA between cells mediated by 
infectious agents such as viruses and conjugative plasmids, respec- 
tively. Because genes involved in these mechanisms reside in the 
genomes of mobile genetic elements, transduction and conjugation 
can be considered side effects of the selfish propagation of these in- 
fectious agents (Redfield 2001). 

Transformation, by contrast, is a bacterium-programmed mech- 
anism of HGT that is mediated by natural competence systems 
encoded in the genomes of many bacteria (Lorenz and Wackernagel 
1994; Claverys et al. 2009). Through transformation, a cell absorbs 
and integrates extracellular DNA (eDNA) into its chromosome. Al- 
though transformation might be viewed as a side effect of metabolic 
consumption of eDNA (Redfield 2001), there is increasing evidence, at 
least in several bacteria, that absorbed eDNA serves as the source of 
information rather than or at least in addition to being a source of 
nutrients (Johnsborg et al 2007; Claverys et al 2009; Zafra et al 2012; 
Johnston et al 2013). Such evidence includes, e.g., the existence of 
competence-induced proteins (e.g., DprA in Streptococcus pneumo- 
niae) that protect absorbed eDNA from degradation by intracellular 
nucleases (for other lines of evidence, see the references cited above). 

Besides transformation, there is another mechanism of HGT called 
gene transfer agents (GTAs), which is encoded in the genomes of 
many prokaryotes (Lang and Beatty 2007; McDaniel et al. 2010; Lang 
et al. 2012). The GTAs are defective virus particles that encapsulate 
apparently random fragments of the host DNA (rather than the genes 
encoding the GTA proteins) and mediate the transfer of these frag- 
ments into other prokaryotic cells. At face value, the GTAs appear to 
function as dedicated vehicles for HGT. The existence of these appar- 
ently prokaryote-programmed mechanisms of HGT suggests an in- 
triguing possibility that these mechanisms have been maintained by 
selection in prokaryotes owing to evolutionary advantages offered by 
HGT — a situation akin to what is currently believed for sexual re- 
combination in eukaryotes (Michod and Levin 1988; Kondrashov 
1993; Keightley and Otto 2006; see Bernstein et al. 1988 and Redfield 
2001 for counteracting views). 

What selective advantages could HGT offer to prokaryotes? To 
address this question, we asked whether an evolutionary advantage of 
recombination that has been originally suggested for eukaryotes could 
be extended to the case of HGT in prokaryotes. In particular, we 
hypothesized that HGT could restore genes that are deleted or 
inactivated by mutations, thereby preventing the stochastic, irrevers- 
ible deterioration of genomes in finite populations known as Muller's 
ratchet (Muller 1964; Felsenstein 1974; Pamilo et al. 1987; Bell 1988; 
Charlesworth et al. 1993; for other types of advantages, see, e.g., 
Redfield et al. 1997; Cooper 2007; Levin and Cornejo 2009; Vos 
2009; Wylie et al. 2010). This hypothesis is consistent with the ob- 
servation of extreme genome reduction and accelerated protein evo- 
lution in intracellular bacterial parasites and symbionts (Moran 1996; 
Novichkov et al. 2009; Merhej and Raoult 2011; McCutcheon and 
Moran 2012). Because of continual bottlenecks and sequestration 
within the host organisms, these bacteria are likely to have small effec- 
tive population sizes and be effectively isolated from HGT between 
populations in different hosts. Under these conditions, the impact of 
Muller's ratchet increases and could lead to genome reduction, accel- 
erated protein evolution, and possibly eventual extinction (Lynch et al. 



1993; Merhej and Raoult 2011). Therefore, HGT might be an important 
requirement for the long-term maintenance of genomic information in 
prokaryotes (Koonin 2011). 

However, HGT in prokaryotes differs from eukaryotic recombi- 
nation both in quantity and in quality. Quantitatively, HGT is 
relatively less frequent and affects a far smaller portion of chromo- 
somes per recombination event than eukaryotic recombination that is 
inextricably linked to sexual reproduction (Dubnau 1999; Feil et al. 
2000). Qualitatively, whereas eukaryotic recombination is a direct ex- 
change of DNA sequences between haploid genomes inherited from 
the two parents, HGT in prokaryotes, in particular, transformation is 
an indirect exchange of DNA sequences between cells via an eDNA 
pool (Redfield 1988). The sources of eDNA can be diverse, but the 
most common one is likely to be the debris of dead conspecific cells. If 
the death rate of a cell increases with the number of deleterious 
mutations in the genome, the genomes of dead cells on average carry 
more deleterious mutations than those of living cells (Redfield et al. 
1997). Consequently, transformation on average is expected to intro- 
duce more deleterious mutations than it would remove. In view of 
these differences between eukaryotic recombination and HGT in pro- 
karyotes, which potentially nullify the hypothetical benefit of HGT, we 
sought to assess the net effect of HGT on the operation of Muller's 
ratchet by performing stochastic population genetic simulations. 

We analyzed a simple evolutionary model of finite prokaryotic 
populations undergoing HGT. Our model is based on the infinite 
population model studied by Redfield (1988) with several extensions 
(see the section Materials and Methods). The model was analyzed 
under the assumption that is most unfavorable to the hypothesized 
advantage of HGT, namely, that variation in fitness is attributed solely 
to variation in death rates, and the fittest individuals do not contribute 
their genomes to the eDNA pool. The results show that, even under 
this assumption, HGT can prevent the operation of Muller's ratchet 
by causing a continual generation of the least-loaded class {i.e., 
genomes with the least number of deleterious mutations) from 
more-loaded classes. Moreover, the prevention of Muller's ratchet 
by HGT is substantially more efficient in a subdivided population 
than in a panmictic population of an equal overall size, provided that 
the diffusion of eDNA is sufficiently rapid. 

MATERIALS AND METHODS 

Definitions for the symbols used in this article can be found in Table 1. 
General description of the model 

The model consists of two components: a population of prokaryotic 
cells and a pool of eDNA. The population undergoes a standard 
mutation-selection process (Wright- Fisher model) as well as HGT, 
through which cells acquire alleles randomly drawn from the eDNA 
pool. The eDNA pool gains input from the population due to natural 
cell death and also undergoes spontaneous decay. The details of each 
of these processes are described in the sections to follow. 

Selection 

The population dynamics was modeled as the discrete-generation 
Wright- Fisher process with fixed population size N. The genome was 
modeled as / loci (haploid), each representing a genomic segment that 
could be replaced by one recombination event (Wylie et al. 2010). The 
typical lengths of such segments are known to be on the order of 10 kb 
in several bacteria (Dubnau 1999; Feil et al. 2000), and the lengths of 
bacterial genomes are on the order of 1 Mb. Accordingly, the value of / 
was set to 100 unless otherwise stated. Each locus was assigned one of 
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Table 1 List of symbols 



Symbol Description 

N Total population size 

U Mutation rate per genome per generation 

s Fitness effect of a mutation (= 0.01) 

r Recombination (HGT) rate per genome per generation 

/ Number of loci in a genome 

deDNA Turnover rate of eDNA per generation 

D p0 p Population migration rate per generation 

D e DNA eDNA diffusion rate per generation 

n. Population size of cells with /' mutations 

m LLC Number of mutations in the least-loaded class 

/rifix Number of mutations fixed in the population 

n mLLC Ne" u / S : the steady-state population size of the least-loaded class 

H mLLC Average gene diversity of the least-loaded class 

H mLLC _i Average gene diversity of one but the least-loaded class 

c/eDNA eDNA quality: the potential of the eDNA pool to remove deleterious mutations in the genomes of the least-loaded class 

teDNA eDNA toxicity: the potential of the eDNA pool to introduce deleterious mutations in the genomes of the least-loaded class 

roi Probability of HGT introducing mutations in the simplified Fisher-Wright model 

r 10 Probability of HGT decreasing mutations in the simplified Fisher-Wright model 

HGT, horizontal gene transfer; eDNA, extracellular DNA . 



the infinite number of possible alleles, each of which was represented 
by a pair of integers: one integer indicating the locus and the other one 
indicating the number of deleterious mutations. All deleterious muta- 
tions were assumed to cause an equal fitness effect s without epistasis; 
thus, the relative fitness f of a genome with i mutations isf = (1 — s) z . 
This simplified model was adopted so that the effect of HGT on 
Muller's ratchet could be assessed without confounding the additional 
effects caused by epistasis (Redfield et al. 1997). The speed of Muller's 
ratchet is known to depend on the two combinations of parameters, 
U/s and sN (Neher and Shraiman 2012). Thus, in this study, the value 
of 5 was set to 0.01, and the values of U and N (more specifically, sN 
exp(— U/s) (denoted as sn mLLC ) and N) were varied to explore the 
behavior of the model. 

Mutation 

In the model, each genome acquired the Poisson distributed number of 
new deleterious mutations with the mean U per generation. Beneficial 
mutations were ignored because we were concerned with the possible 
deterioration or maintenance of initially well-adapted genomes, for 
which the frequency of beneficial mutations is likely to be rare (Maynard 
Smith 1978; Charlesworth et al 1993; cf. Goyal et al 2012). 

Horizontal gene transfer 

HGT generally can be categorized into homologous recombination 
and illegitimate (i.e., nonhomologous) recombination. Because ho- 
mologous recombination is much more frequent than nonhomolo- 
gous recombination (Hulter and Wackernagel 2008; Brigulla and 
Wackernagel 2010), only the former was included in the model. In 
each genome, homologous recombination replaced the Poisson dis- 
tributed number of alleles with mean r per generation per genome 
with alleles randomly drawn from the eDNA pool (after the replace- 
ment the original alleles were discarded). 

Dynamics of the eDNA pool 

The frequency of homologous recombination decreases exponentially with 
the divergence between DNA sequences (Watt et al 1985; Majewski and 
Cohan 1999; Eppley et al 2007). Thus, the most important source of eDNA 
for homologous recombination would be the genomes of conspecific cells. 



For the sake of simplicity, the model assumed that genomes of con- 
specific cells in the same population were the only source of eDNA. 

Prokaryotic cells can release their genomic DNA into the extracellular 
environment either actively by fratricide, suicide, or DNA secretion, or 
passively by natural death. However, the model assumed that genomic 
DNA was released only by natural cell death in order to consider the 
simplest and most unfavorable condition for the putative advantage of 
HGT (see the Introduction section). 

To incorporate the influx of eDNA as the result of natural cell death, 
the model assumed that the per-generation probability of an individual 
with i mutations releasing its genome into the eDNA pool (i.e., dying) 
was Sj = (1 — f ')/ 5^" 0 (1 —fj' )itj where nj is the number of cells with 
j mutations, and fj is the fitness of cells with j mutations normalized by 
the fitness of the least-loaded class, i.e.,fj =fj/f mLLC where m LLC is the 
number of mutations in the least-loaded class (m LLC is defined as the 
integer for which n m = 0 for m < m LLC and n mLLC > 0). 

The aforementioned assumption on 8 Z was obtained from the 
following heuristic argument. Let us consider a continuous population 
dynamics described by dnjd t = bn { — d { n b where b and d { are the birth 
and death rates, respectively. The absolute fitness per unit time of 
a cell with i mutations can be expressed as F { = e h ~ di . The relative 
fitness is expressed as^- = Fj/e b ~ d ° to be consistent with the afore- 
mentioned definition off (so that/ 0 =1). The fraction (probability) of 
cells dying per unit time is 1 — e~ di , which can be expressed as 1 — 
(1 — D 0 )f y where D 0 = 1 — e~ do (D 0 is the probability that an in- 
dividual of the least-loaded class dies). Thus, the normalized contri- 
bution to the eDNA pool by a cell with i mutations is 
(1 — (1 — A) 1/0/2/(1 — (1 ~ Do)fj)rij. Setting D 0 to zero maximizes 
the number of mutations in the genomes of dead cells (see 
Supporting Information, File SI), leading to the expression 
s'j = (1 — //)/2jM)(l ~fj) n j- This expression is used in the seminal 
work of Redfield (Redfield 1988), which considers an infinite popula- 
tion model. In a finite population model, however, the number of 
mutations in the least-loaded class can increase over time owing to 
Muller's ratchet. To keep constant the average number of mutations in 
the eDNA pool relative to the number of mutations in the least-loaded 
class (so that the dynamics is stationary over time), the fitness in s\ 
was normalized by the fitness of the least leaded class f mLLC to obtain 
the expression for 8 Z shown previously. Note that assuming s z is more 
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unfavorable to the hypothesized advantage of HGT than assuming e\ 
because ej is nearly independent of i when m L LC ^> 1- 

In addition, the model assumed spontaneous decay of eDNA at rate 
^eDNA per cell generation. The eDNA decay rate determines the turn- 
over of the eDNA pool, i.e., how synchronous the genetic content of the 
eDNA pool is relative to that of the population (e.g., when ^dna = 1> 
the entire eDNA pool is replaced with the new input every generation). 
Although the absolute concentration of eDNA is likely to influence the 
frequency of HGT (Levin and Cornejo 2009), in the model this in- 
fluence was assumed to be absorbed by r so that the value of *4dna 
does not influence the frequency of HGT. To indicate this fact, ^dna 
will be henceforth referred to as the rate of eDNA turnover. 

With the aforementioned assumptions, the dynamics of the eDNA 
pool was simulated as follows. The state of the eDNA pool was 
described by the vector (e 10 , e n , ■ ■ •eq ' ■ ■ ) where e« is the copy number 
of allele j (i.e.,) mutations) from locus i present in the eDNA pool. In 
each generation, ey reduced by 4dna^ (if 4dna^z)- > 1, it was nor- 
malized to the nearest integer; otherwise, the decrease of e^ was drawn 
from a binomial distribution with the number of trials e^ and success 
probability ^dna)- Then, N cells were randomly chosen with replace- 
ment from the population of the current generation with the proba- 
bility 8f. All alleles in the genomes of the chosen cells were added to 
the eDNA pool, leading to the increase of e y : 

Population structure 

The simplest version of the model assumed that the population is 
genetically isolated from other populations. However, to investigate the 
effect of genetic exchange between semi-isolated populations, we also 
considered an extended model in which the population was subdivided 
into 16 subpopulations of an equal size. The subpopulations were 
arranged in a two-dimensional square grid (4x4 subpopulations) with 
toroidal boundaries. Each subpopulation exchanged individuals with the 
four nearest-neighbor subpopulations via random migration. The num- 
ber of migrants between each pair of neighboring subpopulations was 
drawn from the binomial distribution with the success probability D pop 
and NJ4 trials where N s is the size of a subpopulation (N s = N/16). Each 
subpopulation had its own pool of eDNA that exchanged alleles with 
the four nearest-neighbor pools via diffusion. The copy number of each 
allele diffusing out of a pool was denned as a fraction D eDNA /4 per 
direction per generation (if D eDNA e^/4 > 1, it was normalized to the 
nearest integer; otherwise, the numbers of copies diffusing out were 
drawn from a multinomial distribution with the number of trials e$ 
and success probability D eDNA /4 per direction). 

Organization of the simulation 

The dynamics of the model in one generation was simulated in four 
steps, in the following order: 

1. The Wright- Fisher process. In this step, the population of the current 
generation died and was replaced by that of the next generation. 

2. Production and decay of the eDNA pools. 

3. Diffusion of both the microbial populations and eDNA. 

4. Mutation and HGT. 



RESULTS 

Characterization of Muller's ratchet 

With HGT disabled, the model is known to have the following family 
of stationary distributions «=(«o,«i,"") that are distinguished by 
the number of deleterious mutations in the least-loaded class mux:'- 



h 0 = 0, n mLLC _!=0, n mujC+k =Ne- u / s (U/s) k /k\ where k is 
a non-negative integer (Haigh 1978). A finite-size population cannot 
stay in any of these distributions indefinitely because of stochasticity. 
The absence of back mutations implies that m L Lc increases over time, 
hence the accumulation of deleterious mutations (Haigh 1978). In 
general, such accumulation can result from two different processes 
(Charlesworth et al. 1993): the fixation of deleterious mutations via 
genetic drift or the operation of Muller's ratchet per se, i.e., the extinc- 
tion of the least-loaded class as the result of stochastic fluctuations. 
Muller's ratchet can be reversed by recombination only if a mutation 
has not been fixed (unless eDNA turnover is very slow, as described 
later in this article). Thus, it is beneficial to know through which of 
these processes the accumulation of mutations occurs in the model. 
Thus, we first analyzed the model in the absence of HGT (i.e., r = 0). 
The results of the analysis essentially confirmed the conclusion of 
Charlesworth and Charlesworth (1997) with additional data. Although 
these results in part reproduce the work of Charlesworth and Charlesworth 
(1997), we present them here because they serve as the baseline with which 
to compare the results presented in the later sections. 

The simulations show that the tempo and mode of mutation 
accumulation in the model depend on whether the value of sn mLLC 
(= sNe~ u / s ) is much greater than unity or not, as expected from 
previous work (Rouzine et al. 2008; Neher and Shraiman 2012). For 
sh mLLC = 10 (slow accumulation regime), m LLC and the number of 
mutations fixed in the population (ra^) increase over time in a very 
similar fashion (Figure 1A), so that the number of segregating muta- 
tions in the least-loaded class (m LLC — m^) is at most one. Each 
increase in m LLC is followed by an increase in with a delay of 
approximately a few thousand generations before the next increase 
occurs to m LLC (Figure IB). Thus, the extinction of the least-loaded 
class (i.e., a click of Muller's ratchet) always precedes the fixation of 
a mutation (Charlesworth and Charlesworth 1997). 

The dynamics of mutation accumulation was further analyzed by 
measuring the average gene diversity ("heterozygosity") of the least- 
loaded class (denoted H mLLC ) and that of one but the least-loaded class 
(H muc+1 ; Nei 1987). H miLC is defined as H mLLC = i£| =1 (l - E^x?) 
where Xy is the frequency of allele j in locus i in the least-loaded class 
(Nei 1987). H mLLC+ i is defined in a similar manner. The results show 
that H mLLC and H mLLC+1 display the following dynamics. When Muller's 
ratchet clicks (i.e., when the value of m LLC increases), one but the least- 
loaded class becomes the new least-loaded class by definition, so that 
the value of H mLLC jumps to the value of H mLLC+ i immediately prior to 
the click of the ratchet (Figure 1C). This value of H mLLC+i is always 
nearly 0.02, indicating that, immediately prior to a click of Muller's 
ratchet, one but the least-loaded class consists of / genotypes (one 
mutation in one of the / loci), whose frequencies roughly equal to each 
other (H mLLC+ i = 1 - (1 /if - (l-l//) 2 = 0.0198 for / = 100). Thus, 
immediately after the ratchet clicks, the new least-loaded class consists 
of / genotypes. Because these genotypes have an equal fitness value, the 
gene diversity H mLLC decreases to zero over time owing to genetic drift 
(Figure 1C). Almost simultaneously with H mLLC hitting zero, H mLLC+i 
reaches a steady state value of approximately 0.02 mentioned previously 
(Figure 1C), and the value of increases, indicating the fixation of 
a deleterious mutation (Figure IB). In summary, the accumulation of 
mutations occurs as a cycle of three distinct steps (Figure 2): 

1. Extinction of the least-loaded class due to stochasticity (i.e., Muller's 
ratchet); 

2. Gradual decline of the genetic diversity in the new least-loaded 
class due to genetic drift; 

3. Eventual fixation of a deleterious mutation in the entire population. 
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Figure 1 Operation of Muller's ratchet in the model in the absence of HGT (r = 0). (A, B, D, E) The number of deleterious mutations is plotted 
against time: the population average (gray); the number of mutations fixed in the population (black); and the number of mutations in the least- 
loaded class (red). (C and F) The average gene diversity of the least-loaded class H mLLC (green), that of the one but the least-loaded class H mLLC+ i 
(blue), and the eDNA quality q e DNA (orange) are plotted against time (see the main text for how these quantities are defined). Parameters were as 
follows: s = 0.01, N = 10 5 , d eD NA = 1, / = 100. In the slow ratchet regime, sn mLLC = 10, U = 4.6052 x10" 2 (A-C); in the fast ratchet regime, 
sn mLLC = 1, U = 6.9078X10" 2 (D-F). 



For sn muc = 1, mutations accumulate more rapidly than when 
sn mLLC = 10 (Figure ID). The value of m LLC is incremented multiple 
times before the value of reaches the value of m LLC (Figure IE). This 
result indicates that, unlike the case of sn mi±c = 10, cycles of mutation 
accumulation can overlap with each other in time: a new cycle can begin 
with a click of Muller's ratchet {i.e., step 1 in the previous paragraph) 
before the previous cycle undergoes the fixation of a mutation (i.e., step 
3). Moreover, the dynamics of H mLLC and H miiC+1 are no more syn- 
chronized, in that the decrease of H miLC+1 after a click of Muller's 
ratchet is delayed with respect to the decrease of H miLC (cf. Neher 
and Shraiman 2012). Likewise, the increase of also is delayed with 
respect to the decrease of H miiC+1 . These results indicate that the decline 
of gene diversity after a click of Muller's ratchet (i.e., step 2) is a sequen- 
tial process, in that gene diversity is consecutively lost in different ge- 
notype classes distinguished by the number of mutations {i.e., first, 
H mLLC decreases to zero, and then H mLLC+ i decreases to a steady-state 
level, and so on). Taken together, the results indicate that the 
three- step picture of mutation accumulation described in the 



previous paragraph remains essentially valid in the rapid ratchet 
regime (sh mLLC = 1) as well. 

Therefore, Muller's ratchet is the dominant process leading to the 
accumulation of mutations in the model with HGT disabled (rather 
than the fixation by drift of mutations; Charlesworth and Charlesworth 
1997). However, a click of Muller's ratchet eventually results in the 
fixation of a mutation, so that the number of fixed mutations (m^) far 
exceeds the number of segregating mutations in the least-loaded class 
(^llc — m fhd- This situation is in stark contrast with the result with 
a diploid, sexual model that assumes no recombination and weak 
dominance, which shows that the number of fixed mutations is much 
smaller than the number of segregating mutations (Charlesworth et at. 
1993). The most probable cause of this discrepancy is the prevention of 
mutation fixation in diploid populations due to homozygote disadvan- 
tage, which is absent in our haploid model. The resulting paucity of 
segregating mutations in the haploid model seems to work against the 
putative advantage of HGT; later we show how this limitation might be 
overcome by population subdivision. 



G3'Genes | Genomes | Genetics 



Volume 4 February 2014 I HGT Can Stop Muller's Ratchet I 329 



Least loaded class 



Least loaded class 

r 



Muller's 
ratchet 




Loss of gene 
diversity 



Least loaded class 



Least loaded class 




Figure 2 Schematic diagram depicting accumulation of mutations in 
the model. Horizontal bars indicate the genomes of the least-loaded 
class (genomes of the other classes in the population are not shown). 
Crosses on the bars indicate deleterious mutations. 



Effect of HGT on Muller's ratchet 

We next examined the effect of HGT on Muller's ratchet. For simplicity, 
a rapid turnover of eDNA was assumed by setting 4dna = 1 (the effect 
of slow eDNA turnover is presented in the next section). In this case, 
assuming a very large population, the average number of mutations in 
the eDNA pool per / loci {i.e., per genome) can be approximated by U/s 
+ 1, where U/s is the average number of mutations in the population at 
the steady state with m LLC = 0 (Redfield et al. 1997; File SI). Because the 
mean number of mutations in the eDNA pool is greater than that in the 
population, HGT on average introduces more deleterious mutations 
into the population than it removes and so decreases the steady-state 
population size of the least-loaded class. 

To assess the net effect of HGT on Muller's ratchet, we measured 
the rate of mutation accumulation Am^ I At as a function of 
r (HGT rate) for different combinations of N and sn mLLC values. The 
results show that HGT prevents the accumulation of mutations (i.e., 
Amfix /At = 0) if its frequency is sufficiently high (Figure 3). Moreover, 
the frequency of HGT required to stop Muller's ratchet decreases with 
the increase of N, indicating that the preventive effect of HGT 
increases with N. This trend is seen even in the slow ratchet regime 
(sh mLLC = 10), in which mutation accumulation accelerates with the 
increase of N in the absence of HGT (Figure 3) (note that because 
sn mLLC and 5 are fixed, U was set to a greater value for a greater value of 
N; thus, mutation accumulation can accelerate with the increase of N). 

To understand why HGT can prevent Muller's ratchet despite its 
average deleterious effect due to the high mutation load of eDNA that 
comes from dead cells, we considered the following simple Wright- 
Fisher model. The model assumed two populations: the population of 
the least-loaded class n 0 and that of the other classes n x (the model is 
thus equivalent to one- locus, two-allele, haploid model). The total 
population size N = n 0 + n\ was constant. The value of n 0 at time 
t + 1 was drawn from a binomial distribution with N trials and success 
probability p' denned as follows: 

p' = (1 - 17) (1 - r 01 )p/w + (1 - s')n 0 q/w 

where p = n 0 /N and q = rii/N (= 1 — p) at time t, s' is the fitness 
effect of mutation, U is the mutation rate from n 0 to n x (back 



mutations ignored), w is the average fitness, and r 0 i and r 10 are 
the rates of HGT. For simplicity, HGT was assumed to convert n 0 
into Hi at rate r 01 and, conversely, n x into n 0 at rate r 10 . The effect of 
eDNA was implicitly incorporated into the model by the choice of 
parameters r 0 i and r 10 as follows. If r 0 i and r 10 satisfy the following 
condition, HGT decreases the steady-state level of n Q compared with 
the steady-state level reached in the absence of HGT (i.e., r 0 i = r 10 = 0): 

n 0 (l " s')U'/s' < r 0 i(l - U')(\ - U'/s') 

(the left-hand side of this inequality is the per-generation conversion 
of Hi into n 0 , whereas the right-hand side is that of n 0 into «i, when 
n 0 and n x are set to the steady- state level achieved in the absence of 
HGT). 

With this simplified model, we asked whether HGT facilitated or 
antagonized the survival of n 0 . The mean time to extinction of n Q , 
which is a proxy for Am^/At in the full model, was measured as 
a function of r 0 i and r 10 for different combinations of N and U' values. 
The results show that HGT increases the mean time to extinction of 
n 0 whether or not HGT decreases the steady state level of n 0 , as long 
as the per-generation production of the least-loaded class by HGT 
(ri 0 «i) is on the order of unity or greater (Figure 4; note that n x « N 
for the parameters used in the simulations). This result makes in- 
tuitive sense: the least-loaded class cannot go extinct if it is continually 
produced by HGT in every generation; and as long as such production 
is guaranteed, the steady-state population size of n 0 is immaterial to 
the extinction of n 0 . To check whether this explanation also applies to 
the full model, we investigated the consistency between the two mod- 
els from a few different angles as described below. 

First, the simplified model indicates that the required frequency of 
HGT to prevent Muller's ratchet decreases as N increases because 
increasing N increases the production rate of the least-loaded class 
(with ho kept constant). The same result is produced by the full model 
as already described previously. 

Next, Figure 4 shows that whether HGT increases the time to 
extinction depends on the value of r 10 much more strongly than on 
the value of r 01 (as seen from the fact that the boundary between 
yellow and red regions is nearly horizontal). We thus hypothesized 
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Figure 3 The effect of HGT on the operation of Muller's ratchet. The rate 
of mutation accumulation Amf ix /At is plotted as a function of the rate of 
HGT r. The accumulation rate was estimated as m^-Jt, where mf ix is the 
number of fixed mutations at the end of a simulation and t is the duration of 
a simulation (e.g., t > 1 0 6 for N = 1 0 6 and t > 1 0 9 for N = 2 x 1 0 3 ). The error 
bars indicate standard error of means estimated as ^/m^/t 2 . The param- 
eters were as follows: s = 0.01, d e DNA = 1, / = 100. In the slow ratchet 
regime, sn mLLC = 10 (filled circles): N = 10 6 and U = 6.9078 x 10" 2 (black), 
N = 10 5 and U = 4.6052 x 10" 2 (red), N = 10 4 and U = 2.3026 x 10" 2 
(green), N = 2 x 10 3 U = 6.931 5 x10" 3 (blue). In the fast ratchet regime, 
sn mLLC = 1 (open circles): N = 10 6 and U = 9.2103x10" 2 (black), N = 10 5 
and U = 6.9078 x10" 2 (red), N = 10 4 and U = 3.6052 x 10" 2 (green), 
N = 2 x 1 0 3 and U = 2.9957 x 1 0" 2 (blue). 



that, in the full model, HGT prevents Muller's ratchet most effectively 
when the probability of HGT producing the least-loaded class is max- 
imized. This probability can be approximated by the probability P 10 
that an HGT event transforms one but the least-loaded class 
(n mLLC - i) into the least-loaded class (n mLLC ), because the contributions 
from other mutant classes likely involve more than one HGT event. 
P 10 can be calculated as — l~ 1 )e~ l u l s under the assumption 
that N is very large (File SI). P 10 takes the maximum value when 

r l = P 1 -5/1/ +1/2- ^(s/U) 2 + 1/4 (F l ks/U + (s/U) 2 

when s/U <C 1; File SI). The non-monotonic dependency of P 10 on 
/ can be intuitively understood as follows: decreasing / can decrease P 10 
because it increases the probability that an allele randomly drawn 
from the eDNA pool contains a mutation that is not carried by the 
least-loaded class (if / = 1, this probability is one). However, increasing 
/ also can decrease P 10 because it decreases the probability that HGT 
occurs to the very locus in which one but the least-loaded class has 
the mutation not carried by the least-loaded class. Therefore, P 10 
takes the maximum value at an intermediate value of /. By contrast, 
the probability P 01 that HGT introduces one or more mutations 
in the least-loaded class is 1 — (1 — l~ 1 )e~ l u / s (File SI), which is 



a monotonically decreasing function of / because increasing / decreases 
the probability that an allele randomly drawn from the eDNA pool 
contains mutations. Thus, the aforementioned hypothesis can be 
reformulated as follows: Amfix/At reaches the minimum value when 
/ assumes the value/. To examine this hypothesis, we measured Amfe/ 
At as a function of / for various combinations of N and sn 0 (Figure 5). 
The results show that Am^/At takes the minimum value when the 
value of / is on the same order of magnitude as / (around 10 rather 
than 100), in accord with the hypothesis. 

Finally, the simplified model shows that r^ni > 1 is sufficient for 
HGT to prevent the extinction of no- To check if a similar condition 
applies to the full model, we calculated the per-generation production 
of the least-loaded class by HGT occurring to one but the least-loaded 
class (rPio«i). For sn mLLC = 10, Am^c/ At sharply drops at r « 0.01 
for N = 10 6 and 10 5 (Figure 3). The values of rPi 0 n mLLC in these cases 
are approximately 0.6 and 0.4 for N = 10 6 and 10 5 , respectively, in 
qualitative agreement with the expectation that rPiori mLLC be on the 
order of unity (note that rPi 0 n mLLC < 1 is expected because the con- 
tributions from the more-loaded classes are ignored and n mLLC > 0 is 
a stronger condition than Arafi x / At = 0). For sn mLLC = 1, Am LLC /At 
sharply drops at r « 0.1 for N = 10 6 and 10 5 . The values of rPi 0 n mLLC 
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Figure 4 The mean time to extinction of the least-loaded class in the simplified Fisher-Wright model. Below the diagonal lines, HGT decreases 
the steady-state population size of the least-loaded class (i.e., the condition rio(1 - s')U /s' < roi (1 - l/')(1 - U'/s') is fulfilled). The mean time to 
extinction is normalized by that obtained in the absence of HGT (i.e., roi = no = 0). The parameters were as follows: s' = 0.01, and s'no = 10 
where n 0 = N(1 - U'/s'). N= 10 4 and U' =9x10" 3 (A); N= 10 5 and U' =9.9x10" 3 (B); N= 10 6 and U' = 9.99x10" 3 (C). Each data point was 
obtained as an average of 100 simulation runs. 
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Figure 5 The normalized rate of mutation accumulation as a function of 
the number of loci /. The reference rate is set to the rate obtained for / = 
100. The value of / that maximizes the probability P 10 (that HGT con- 
verts one but the least-loaded class to the least-loaded class) is indi- 
cated as / in the graph. The parameters and color coding are the same 
as in Figure 3 (except for /)• 
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in these cases are approximately 0.8 and 0.6 for N = 10 6 or 10 5 , 
respectively, again in qualitative agreement with the expectation. 
(The results were also consistent for smaller values of N, at least in 
the orders of magnitude.) 

Taken together, the consistency of the results indicates that the 
simplified model captures, at least qualitatively, the mechanism by 
which HGT prevents Muller's ratchet in the full model. Namely, HGT 
can prevent Muller's ratchet even if it decreases the steady-state pop- 
ulation size of the least-loaded class because it can hinder the extinc- 
tion of the least-loaded class by continually generating it from the 
more-loaded classes. 

Effect of the eDNA turnover rate on the prevention of 
Muller's ratchet by HGT 

We further investigated the effect of the eDNA turnover rate on the 
ability of HGT to prevent Muller's ratchet. As described previously, 
the accumulation of mutations in the model is initiated by a click of 
Muller's ratchet followed by the fixation of a mutation. Once a muta- 
tion is fixed (i.e., once a good allele is lost), HGT cannot undo it 
because the only source of HGT in the model is eDNA derived from 
the members of the same population. However, slow turnover of 
eDNA can delay the loss of good alleles from the eDNA pool and 
thereby might help HGT prevent the accumulation of mutations. 

To examine this possibility, we measured Am^/ At as a function 
of 4dna for various combinations of N and sn mLLC values (Figure 6). 
The results show that Am^/At decreases with the decrease of rf eDNA . 
However, for this effect to be significant, eDNA turnover has to be 
slower than the population turnover (i.e., generation time) by more 
than three orders of magnitude (4dna < 10~ 3 ), a condition that 
appears to be unrealistic in natural environments (Nielsen et al. 
2007; Pietramellara et al. 2009). 

To understand why such a severe condition applied, we 
measured the "quality" of eDNA denned as ^dna = 
T!k=iT,jLoY,i°=j+i( i ~ j) x kiykj where x ki is the fraction of the 
least-loaded class having allele i (i.e., i mutations) in locus k, and 
ykjis the fraction of allele j from locus k in the eDNA pool. g e DNA 
indicates the potential of the eDNA pool to remove deleterious 
mutations in the genomes of the least-loaded class (e.g., g eDNA is 
zero if the eDNA pool contains no allele that can remove mutations 
in the least-loaded class). g eDNA was measured in the absence of 
HGT for various values of 4dna- The results show that, in the slow 



ratchet regime (sh 0 = 10), the dynamics of ^ eDNA closely follows 
that of H mLLC when eDNA turnover is not very slow (4dna — 
10~ 2 ) (Figure 1C and Figure 7, A and B). When eDNA turnover 
is very slow (4dna — 10~ 3 ), the decay of ^dna lags behind the 
decline of H mLLC (Figure 7C). This result indicates that eDNA turn- 
over has to be slower than the decline of H mLLC to delay the decay of 
^eDNA significantly. Because the decrease of H mLLC is attributed to 
genetic drift as mentioned previously, its timescale is determined 
by the population size of the least-loaded class, which was approx- 
imately 1000 in the simulations for the slow ratchet regime 
(sh mLLC = 10 and 

n m LLC — 1000)- Thus, 4dna must be <10 3 to 
cause any significant effect. 

The aforementioned argument might imply that the effect of slow 
eDNA turnover could be more significant if the population size of the 
least-loaded class is small. However, this turns out not to be the case. 
In the fast ratchet regime (sn mLLC = 1 and h mLLC = 100), the dynamics 
of ^eDNA no more closely follows the dynamics of H mLLC> and it still 
takes more than 500 generations for g eDNA to decrease to zero when 
eDNA turnover is rapid (Figure 7, D— F). Therefore, 4dna still has to 
be smaller than 1/500 to have any significant effect (File SI contains 
a more detailed explanation of the results described in this and the 
previous paragraphs). 

In summary, the fixation of mutations is an evolutionary process 
that is driven by selection and drift and thus occurs on relatively long 
timescales compared with the generation time of individuals, so that 
eDNA turnover has to be commensurately slow to have any 
appreciable effect. 

Effect of population subdivision on the prevention of 
Muller's ratchet by HGT 

Next, we considered population subdivision. On the one hand, 
population subdivision can accelerate accumulation of mutations 
because it enhances the effect of genetic drift as shown in previous 
studies (Higgins and Lynch 2001; Combadao et al. 2007). On the 
other hand, subdivision might help HGT prevent Muller's ratchet as 
follows. Let us suppose that the isolation between subpopulations is so 
strong that the accumulation of mutations in different subpopulations 
is independent of each other in the absence of HGT. In this case, 
a mutation must be fixed in every subpopulation independently before 
being fixed in the entire population. This situation is likely to lead to 
an increased number of segregating mutations when the entire 
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Figure 6 Effect of slow eDNA turnover on the prevention of Muller's 
ratchet by HGT. The rate of mutation accumulation Amf ix /At is plotted 
as a function of the eDNA turnover rate d eDNA . The parameters and 
color coding are the same as in Figure 3 (except for d eDNA ). 



population is considered (rather than when each subpopulation is 
considered separately). The increased number of segregating muta- 
tions by itself does not affect accumulation of mutations if recombi- 
nation requires direct contact between individuals as in many animals. 
HGT, however, does not require such direct contact. If subpopulations 
share a common eDNA pool because of the rapid transport of DNA 
molecules between their habitats, increasing the number of segregat- 
ing mutations is expected to increase the number of mutations that 
are introduced or removed by HGT events. This should enhance the 
preventive effect of HGT on Muller's ratchet because removal of 
mutations by HGT is more important than introduction of mutations 
by HGT according to the simplified Wright- Fisher model presented 
above. The key question is whether this potential advantage of pop- 
ulation subdivision can outweigh the disadvantage of population sub- 
division associated with the enhanced genetic drift. 

To address the aforementioned question, we measured Amfix/At 
as a function of population migration rate (D pop ) in a model incor- 
porating population subdivision (Materials and Methods), was 
denned as the average of measured within each subpopulation. 
For simplicity, all subpopulations were assumed to share a common 
eDNA pool (i.e., D eDNA = oo ; we will consider finite D eDNA later). The 
results show that, in the absence of HGT (r = 0), the system becomes 
more prone to Muller's ratchet (i.e., Amfi x /At increases) as D pop 
decreases (Figure 8, A and B). This outcome is expected, given the 
disadvantage of population subdivision due to enhanced genetic drift. 
By contrast, in the presence of HGT (r > 10~ 4 ), the dependency of 
Amfix/At on D pop is non-monotonic: the system initially becomes 
more prone to Muller's ratchet as D pop decreases, but becomes less 
so as D pop decreases further, with the position of the maximum 
Amfix/At along the D pop axis shifting as a function of r (Figure 8, 
A and B). The comparison against the curve for r = 0 shows that the 
reduction in Amfix/At caused by HGT increases as D pop decreases, 
indicating that population subdivision increases the effect of HGT. At 
a high HGT rate (r = 0.01), in the fast ratchet regime (sh mLLC = 1), 
the system is more resistant to Muller's ratchet at lower values of 
Dp 0 p (<10~ 5 ) than at greater values of D pop (>10 _1 ) for which the 
system can be considered well-mixed (Figure 8B). In the slow ratchet 
regime (sh mLLC = 10), however, this situation cannot be observed be- 
cause Amfix/At falls to zero and so shows no difference between 
low and high values of D pop within the accuracy of the simulations 
(Figure 8A). In addition to the aforementioned results, Amfi x /At was 
measured as a function of D pop for other combinations of N and 



sn m LLC values; the results were qualitatively the same as described 
previously (File SI). Taken together, these results indicate that the 
advantage of population subdivision can outweigh the associated dis- 
advantage, provided that the frequency of HGT is sufficiently high, 
and subpopulations share a common eDNA pool. 

To elucidate the mechanism by which population subdivision 
helps HGT prevent Muller's ratchet, we first sought to test the expec- 
tation that population subdivision increases the number of segregating 
mutations at the level of the entire population. To this end, we mea- 
sured the average gene diversity of the least-loaded class within and 
between subpopulations (Nei 1987). The average gene diversity within 
one subpopulation is denned as Hs y = (l/OSjLif^ — ^2^=o x lki) 
where x^i is the frequency of allele i in locus k in the least-loaded 
class of subpopulation y (the least-loaded class was denned separately 
for each subpopulation). The average gene diversity within subpopu- 
lations H s is denned as the average H Sy over all subpopulations 
(weighted by the population size of the least-loaded class in each 
subpopulation). The average gene diversity in the entire population 
H T is denned similarly with x y ^ replaced by the frequency of an allele 
in the least-loaded classes from all subpopulations combined together. 
The average gene diversity between subpopulations was denned as 
H T — H s (Nei 1987). H T and H s were measured in the absence of 
HGT. Simulations show that the values of H T and H s initially increase 
and then saturate over time (results not shown). The time averages of 
H T and H s after or near saturation were obtained as a function of D pop 
(Figure 9A). The result shows that H T increases as D pop decreases, 
whereas H s remains close to zero for the entire range of D pop . This 
result indicates that population subdivision increases the number of 
loci that are polymorphic between subpopulations, whereas within 
subpopulations almost all loci remain monomorphic. Therefore, pop- 
ulation subdivision increases the number of segregating mutations at 
the level of the entire population as expected. 

The increase in the number of segregating mutations is expected to 
drive an increase in the number of mutations that are introduced or 
removed by HGT events. To test this expectation, we measured the 
quality of the eDNA pool denned above (^dna) and the "toxicity" of 
the eDNA pool denned as f e DNA = SLiS^iSm^O' ~ 0*^; 
(<2eDNA and t eD NA were calculated for each subpopulation and aver- 
aged over subpopulations using the population sizes of the least- 
loaded classes as weights). t eDNA indicates the potential of the eDNA 
pool to introduce deleterious mutations in the genomes of the least- 
loaded class. Note that the difference t eDNA — g eDNA is equal to the 
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Figure 7 Effect of slow eDNA turnover on the quality of eDNA q e DNA • The average gene diversity of the least-loaded class H mLLC (green), that of 
the one but the least-loaded class H mLLC+ i (blue), and the quality of eDNA q e DNA (orange) are plotted against time for various values of eDNA 
turnover rate d eDNA . The parameters were as follows: / = 100, s = 0.01, c/ s dna = 0-1 (A and C), 0.01 (B and E), and 0.001 (C and F). In the slow 
ratchet regime, sNe'^ 5 = 10, N = 10 5 and U = 4.6052 x 10" 2 (A-C). In the fast ratchet regime, sNe" u / s = 1 , N = 1 0 5 , and U = 6.9078 x10" 2 
(D-F). 



average number of mutations introduced (or removed if the value is 
negative) by one HGT event occurring to an individual of the least- 
loaded classes. Thus, — g eDNA can be interpreted as the contribution 
to the average number of mutations introduced by HGT by those 
HGT events that decrease the number of mutations, whereas ( eDNA 
as the contribution by those events that increase the number of muta- 
tions. Simulations show that the values of <jeDNA and ( c dna increase 
and then saturate over time except for D pop = 0 (result not shown; for 
D pop = 0, saturation is likely to be reached, but probably requires an 
exceedingly long time). The time averages of <2eDNA and ^dna after 
saturation were obtained as a function of D pop in the absence of HGT 
(Figure 9B; for D pop = 0, the values reached toward the end of the 



simulation are shown, thus indicating a lower bound). The result 
shows that population subdivision increases the number of mutations 
removed by an HGT event (^ eDNA ) as well as the number of muta- 
tions introduced by an HGT event (^dna)> thus confirming the ex- 
pectation (Figure 9B). Moreover, it leaves nearly constant the average 
(i.e., net) number of mutations introduced by an HGT event (f eDNA — 
g eDNA ). In other words, population subdivision increases the conver- 
sion of the least-loaded class by HGT in the direction of both in- 
creasing and decreasing the number of mutations without changing 
the magnitude and direction of the net conversion. According to the 
simplified Wright-Fisher model described before, removal of muta- 
tions by HGT (r 10 ) has a greater influence on the extinction of the 
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Figure 8 Effect of population subdivision on the prevention of Muller's ratchet by HGT when eDNA diffusion is infinitely fast (D eDNA = qo). The 
rate of mutation accumulation Arhf jx /At is plotted as a function of population migration rate D pop for various values of HGT rate r. The color coding 
is as follows: r = 0 (black), r = 10~ 4 (red), r = 10~ 3 (green), and r = 10~ 2 (blue). The parameters were as follows: in the slow ratchet regime, 
s N e -u/s = 10, N= 10 5 , and U = 4.6052 x 1 0" 2 (A) ; in the fast ratchet regime, sNe" u / s = 1, N= 10 5 , and U = 6.9078 x 1CT 2 (B). The number of 
subpopulations was 4x4 with toroidal boundaries (see Materials and Methods). 



least-loaded class than introduction of mutations by HGT (r 0 i). Al- 
though r 10 and r 0i are concerned with HGT occurring in different 
genotypes, they are analogous to ^dna and ( eDNA , respectively. This 
analogy points to the advantage of increasing g eDNA and f eDNA with 
^dna — ^eDNA kept constant to prevent Muller's ratchet. 



Whether population subdivision is advantageous or disadvanta- 
geous depends on the combination of parameters. To examine what 
appears to be the most important aspect of this parameter de- 
pendence, we next considered the finite diffusion of eDNA. To this 
end, the value of Amixc/Af was compared between D pop = 10~ 5 
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Figure 9 Effect of population subdivision on the genetic structure of 
populations and on the contents of eDNA pools. (A) The average gene 
diversity within subpopulations (H s ) and in the entire population (H T ), 
and the coefficient of gene differentiation G S t = (Hj - H S )/H T are 
plotted as a function of D pop in the absence of HGT (r = 0). The plotted 
values were obtained as the time average after or nearly after the values 
reached saturation (all time averages were taken from the last 2x10 4 
generations of simulations). The error bars denote standard deviation. 
The parameters were the same as in Figure 8A (filled circles) and Figure 
8B (open circles) except that r = 0. (B) The quality q e DNA and toxicity 
teDNA of eDNA are plotted as a function of D pop in the absence of HGT 
(r = 0). The values of q e DNA and t e DNA were obtained as the time 
average after the values reached saturation except for D pop = 0, for 
which the values indicate lower bounds (all time averages were taken 
from the last 2 x 1 0 5 generations of simulations). The error bars denote 
standard deviation. The parameters were the same as in Figure 8B 
except that r = 0. 
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(subdivided population) and D pop = 10 -1 (well-mixed population) 
for various combinations of 4dna and D eDNA values. As expected, 
the results show that Amixc/Af is smaller for D pop = 10~ 5 than for 
D pop = 10 _1 when <i e DNA is sufficiently small or D eDNA is sufficiently 
large (Figure 10). The values of D eDNA and 4dna demarcating the 
region of the parameter space for which population subdivision is 
advantageous do not appear unrealistic, although it is difficult to judge 
whether the conditions represented by this region are achievable in 
natural environments (see also Discussion). 

DISCUSSION 

The results presented here suggest that HGT can prevent the 
operation of Muller's ratchet in prokaryotic populations, even if on 
average HGT introduces more deleterious mutations than it removes. 
The avoidance of Muller's ratchet via transformation and recombina- 
tion with eDNA might appear somewhat paradoxical because on 
average eDNA is expected to carry more deleterious mutations than 
the DNA in live prokaryotic cells. Indeed, it has been long recognized 
that "sex with dead cells" is a dubious proposition (Redfield 1988, 
2001; Redfield et al. 1997). However, the modeling results indicate 
that, even though transformation via eDNA is expected to increase the 
mean mutation load and hence decrease the average fitness of a pop- 
ulation, it nevertheless can stop Muller's ratchet. This appears to be 
the case because HGT provides for the chance to eliminate deleterious 
mutations, leading to the continual restoration of the least-loaded 
class. In other words, HGT prolongs the persistence of the least- 
loaded class in the population in the face of stochastic fluctuations, 
hence the prevention of Muller's ratchet. 

The efficacy of HGT in the prevention of Muller's ratchet depends 
on the stability of eDNA, its diffusion rate, and population subdivi- 
sion. The contribution of eDNA stability (i.e., slow turnover) is in- 
tuitively clear because it delays elimination of high-fitness alleles from 
the population and so increases the chance that these alleles offset 
Muller's ratchet via HGT. However, the simulations show that, for this 
effect to be substantial, the characteristic lifetime of eDNA should 
exceed the generation time by orders of magnitude (Figure 6), sug- 
gesting that eDNA stability alone is not a key factor in the evolution of 
microbial populations. If, however, a population is subdivided, en- 



hanced eDNA stability facilitates the transport of eDNA across sub- 
populations and thus contributes to the ability of HGT to prevent 
Muller's ratchet. Under these conditions, for the effect to be signifi- 
cant, the characteristic lifetime of eDNA does not have to be exceed- 
ingly long, depending on the diffusion rate of eDNA (Figure 10). 

If the transport of eDNA is sufficiently fast, HGT can prevent 
Muller's ratchet more efficiently in a subdivided population than in an 
undivided population. As a result, a subdivided population can be 
more resistant to Muller's ratchet than an undivided population of 
an equal overall size despite the fact that population subdivision accel- 
erates Muller's ratchet in the absence of HGT owing to the enhanced 
effect of genetic drift. Put differently, to maintain genomic informa- 
tion in the face of deleterious mutations and genetic drift, it is advan- 
tageous to partition individuals into multiple (smaller) subpopulations 
and let them "cross-reference" each other's genetic information 
through recombination rather than collect all individuals in one pop- 
ulation and thereby maximize the efficacy of natural selection. Note, 
however, that although population subdivision can be beneficial, it is 
not necessary for HGT to prevent Muller's ratchet. 

How do these findings relate to the previous population genetics 
studies on the advantage of recombination? Previous studies have 
shown that recombination is beneficial (compared with no recombi- 
nation) when a population suffers from the reduced efficacy of natural 
selection because of negative linkage disequilibrium (i.e., biased asso- 
ciations on chromosomes between beneficial and deleterious alleles) 
which arises from genetic drift (the Hill-Robertson effect) or syner- 
getic epistasis (Hill and Robertson 1966; Felsenstein 1974; Feldman 
et al. 1980; Kondrashov 1988, 1993; Otto and Lenormand 2002). Our 
present results are consistent with these findings because they show 
that the reduction in the rate of mutation accumulation caused by 
HGT is greater in a subdivided population than in an undivided 
population (Figure 8, A and B). Indeed, population subdivision causes 
negative linkage disequilibrium in many pairs of loci in the least- 
loaded classes at the entire population level because the least-loaded 
classes in different subpopulations have mutations in different loci 
(this can be seen from the fact that H T > 0 and H s « 0 when 
^pop ~ 0)- I n addition, the results show that the rate of mutation 
accumulation itself can be lower in a subdivided population than in an 
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Figure 10 Effect of population subdivision on the prevention of 
Muller's ratchet by HGT when the eDNA diffusion rate is finite. Red 
squares indicate the combinations of c/ s dna and D eDNA values for which 
the rate of mutation accumulation Arhf ix /At is lower when a population 
is subdivided (D pop = 10~ 5 ) than when a population is well-mixed 
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10 1 ). Black circles indicate the parameter region where the 



opposite was the case. The parameters (other than D eDNA and 2 eDNA ) 
were the same as in Figure 8B with r = 0.01 . 
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undivided population, despite the disadvantage of population subdi- 
vision due to enhanced drift. This conclusion does not immediately 
follow from the previous findings because the fact that negative link- 
age disequilibrium increases the benefit of recombination (compared 
with no recombination) does not necessarily mean that negative link- 
age disequilibrium is beneficial in preventing the accumulation of 
mutations in the presence of recombination. 

Evolutionary maintenance of genomic information based on 
population subdivision and recombination hinges on the fact that 
genetic exchange between cells happens via a shared pool of eDNA. 
On the one hand, this implies that the indirectness of HGT might 
actually be an asset rather than a drawback for prokaryotes because 
HGT provides the possibility of exchange between physically separated 
subpopulations that accumulate different mutations. 

On the other hand, indirect recombination imposes strict 
conditions under which population subdivision is advantageous in 
terms of Muller's ratchet, requiring restricted migration between pop- 
ulations and rapid transport of eDNA molecules. Although at present 
it is difficult to judge whether such conditions are fulfilled in natural 
environments, the available data seem to suggest soils as an environ- 
ment that is conducive to these processes. The distribution of bacteria 
in natural soils is "patchy," i.e., bacteria typically occur as small col- 
onies separated from each other, and their locations are restricted by 
the availability of soil micro-pores, water, and organic substances 
(Foster 1988; Stotzky 1989; Grundmann 2004; Young et al. 2008; 
Vos et al. 2013). Moreover, bacterial motility can be highly restricted 
by water availability (Dechesne et al. 2010). In addition, soil bacteria 
display high levels of genetic diversity within local populations 
(Grundmann 2004). For example, Agrobacterium biovar 1 and Nitro- 
bacter-like bacteria within minute soil samples (<1 mm in diameter) 
display genetic divergence as great as that between strains sampled 
from different geographical areas (Grundmann and Normand 2000; 
Vogel et al. 2003, although Agrobacterium biovar 1 is not known to be 
naturally competent, this bacterium readily undergoes homologous 
recombination; Costechareyre et al. 2009). Also, genomes of Pseudo- 
monas stutzeri display extremely high gene diversity (H T > 0.8) and 
high frequency of null alleles (>80% of the examined strains fail to 
exhibit the activity of at least one of the 20 enzymes examined; Rius 
et al. 2001). Taken together, these findings suggest that bacterial pop- 
ulations in unsaturated soils are strongly subdivided and genetically 
heterogeneous within local environments. This population structure is 
compatible with a situation in which HGT and population subdivision 
jointly prevent the operation of Muller's ratchet (if eDNA transport is 
sufficiently rapid as well). 

Moreover, eDNA in soils becomes highly stable against hydrolysis 
by nucleases when it is absorbed to clay minerals, sand particles, and 
humic acids without losing the ability to transform competent bacteria 
(Vries and Wackernagel 2005; Nielsen et al. 2007; Pietramellara et al. 
2009). Natural transformation assays indicate that DNA added to 
nonsterile soils retains transformation ability for 3—15 days (Gallori 
et al. 1994; Sikorski et al. 1998). In addition, eDNA can be transported 
within unsaturated soils through water capillarity or leaching (Ceccherini 
et al. 2007; Pote et al. 2010). Although the difference between bacteria 
and eDNA in terms of transport efficiency within soils remains to be 
evaluated, the obvious size difference suggests that eDNA is more readily 
transported than bacteria. Taken together, these findings suggest soil 
prokaryote systems as a prime candidate for a system in which HGT 
and population subdivision could be important for the long-term main- 
tenance of genomic information in the face of Muller's ratchet. Further- 
more, it appears plausible that both the molecular machinery of natural 
transformation and GTAs evolved, at least in part, as adaptations 



facilitating HGT and overcoming the limitations on HGT imposed 
by the aforementioned conditions, e.g., GTAs can increase the stability 
of DNA emitted from the cell. 

Is Muller's ratchet an important factor of evolution in prokaryotes? 
Its role can be questioned because many bacteria have large popula- 
tion sizes in which the extent of stochasticity is limited and that are 
accordingly immune to Muller's ratchet (but see also the discussion of 
soil bacteria above in this section). However, bacteria are also known 
to go through frequent population bottlenecks under stress conditions 
(Levin et al. 2000; Aertsen and Michiels 2004; Wolf et al. 2005), and it 
is under these conditions that Muller's ratchet and mechanisms that 
can overcome it could become relevant. Indeed, competence has been 
shown to increase stress resistance in bacteria (Engelmoer and Rozen 
2011) and induction of the competence regulons is an important part 
of bacterial stress response (Claverys et al. 2006; Prudhomme et al. 
2006). Although the biology of the GTAs is not nearly as well un- 
derstood as that of competence and transformation, stress -induced 
prophage replication that is likely to potentiate HGT has been dem- 
onstrated (Garcia- Russell et al. 2009). 

To conclude, the results of the modeling study described here 
indicate that HGT can prevent the operation of Muller's ratchet, and 
this effect can be enhanced by population subdivision. Thus, through 
preventing the action of Muller's ratchet, HGT could be a necessary 
condition for the long-term survival of prokaryotic populations. The 
benefits of HGT for prokaryotes certainly are not limited to the pre- 
vention of Muller's ratchet and additionally include acquisition of 
novel genes and functions (see, e.g., Levin and Cornejo 2009; Wylie 
et al. 2010). Jointly, these advantages could be the driving forces 
behind the evolution of dedicated mechanisms for HGT. 
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